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Motivation 


Introduction and Motivation 
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Where is the finite-density QCD? 


Neutron star 

e Cold and dense nuclear matters 

o 2Msun neutron star (2010) 

o Gravitational-wave observations (2017~) 
Heavy-ion collision 


e Low-energy scan of heavy-ion collisions are planned and run in 
many facilities (RHIC, SIS, J-PARC) 


Reliable theoretical approach to equation of state must be developed! 


A(T, u) = fra Det (P(A) + m) exp —Sym(A). 
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Silver Blaze problem in finite-density QCD 
QCD & |QCD| 


ZQCD = [PA (dety D, ) e 9v. Z|QCD| = [ra Idet^, D, | e ?YM. 


At u = 0, these two are the same! But, 


ng 
QCD 


1QCD] 


m /2 my/3 "H 


At T = 0, the state must be equal to the QCD vac. for |u| < my/3. 
Can we describe this "trivial" behavior using path integral? 
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Lefschetz-thimble integrations 


Lefschetz decomposition formula 


Oscillatory integrals with many variables can be evaluated using the 


“steepest descent” cycles J,: (classical eom S'(2,) = 0) 


f d"z e 3) — Y (Ke; y) d'zg 9i 


o 


Jz are called Lefschetz thimbles, and Im[S] is constant on it: 


M =}, dzi(t) _ eo 


1—-—oo dt pe 


= [40 


(Ko, R): intersection numbers of duals K, and R” 


(Ko = {2(0)|2(00) = 24 )). 


[Witten, arXiv:1001.2933, 1009.6032] 
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Lefschetz-thimble integrations 


Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) = | evi (5 + 2 


Complexify the integration variable: z = x + iy. 


Im(¢) 


4 


Integrand on R, and on A 
(a — 1) 
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Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (c = 1,2) for the Airy 


integral: 
d 3 
Ai(a) = > No [ Pi (5 taz) : 


ng: intersection number of the steepest ascent contour Ke and R. 


Figure: Lefschetz thimbles J and duals K (a = 100%, =i) 
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Sign problem & One-site Hubbard model 


Simplest model study to understand the sign problem 
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Sign problem & One-site Hubbard model 


One-site Fermi Hubbard model 
One-site Hubbard model: 
A = Ufuü, — u(fu + fy). 
Fock state gives the number density immediately: 


10 2(e9" + ePQu-U)) 
bôu id= 1 + 2e8Hn + ef(2u-U)" 


In the zero-temperature limit, 
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Sign problem & One-site Hubbard model 


Path integral for one-site model 


The path-integral expression for the one-site Hubbard model: : 


B n i 2- UD 
a=a > | ay (14 ef Ge)? 78/20 
27U R 2 Ub cia addio D 


Fermion Det 


Integrand has complex phases causing the sign problem. 


( is an auxiliary field for the number density: 


Yuya Tanizaki (University of Tokyo, RIKEN) Lefschetz-thimble path integral Nov 25, 2015 @ KEK 


10/22 


Sign problem & One-site Hubbard model 


Flows at y/U < —0.5 and u/U > 1/5 


Re(z) 


Figure: Flow at u/U 


2 


| dz e 8), 
Se 


Number density: n, = 0 for y/U < —0.5, n, = 2 for u/U > 1.5. 


(YT, Hidaka, Hayata, 1509.07146) 
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Lefschetz-thimble path integral 
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Sign problem & One-site Hubbard model 


Flows at —0.5 < u/U < 1.5 
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Figure: Flow at y/U = 0 


py p 1. 


This value is far away from n = Im (z)/U=0, 1, or 2. 


ra 


Complex saddle points lie on Im(z,,)/U 


12722 


Nov 25, 2015 @ KEK 


Lefschetz-thimble path integral 
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Curious incident of n in the one-site model 


We have a big difference between the exact result and naive 
expectation: 


05 0 0.5 1 15 uU 


(YT, Hidaka, Hayata, arXiv:1509.07146[hep-th]) 
(cf. Monte Carlo with one-thimble approx. gives the naively expected results. 


Fujii, Kamata, Kikukawa,1509.08176, 1509.09141; Alexandru, Basar, Bedaque,1510.03258) 


Question 
Can we unveil this curious incident? 
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Semiclassical analysis 


Semiclassical study of the curious incident 
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Semiclassical analysis 


Complex classical solutions 
Let us solve, for —0.5 < u/U < 1.5, 


B. 2B exp B (iz + +4) 
iz 
U 1 - exp 8 (iz - uu 4- E) 


If GU > 1, the solutions are labeled by m € Z: 
U 
Zm 221 (; + 3) + 27 m1. 


At these solutions, the classical actions become 


ON 


2 2 
Re (Sm — So) = zm, 
Im Sm c 27m (5 4 1). 
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Semiclassical analysis 


Semiclassical partition function 
Using complex classical solutions Zm, let us calculate 


oo 
La := > p nm. 


Tn—-—oo 


This expression is valid for —1/2 S u/U $ 3/2. 
This is calculable using the elliptic theta function: 


Za 


K 


= 1 n 
-S [149 2 Heo +1 ntm?/pu 
e + 2 Tm (5 + 2 e 


— eo pe 1 —2n? /BU 
e o ((5*5).« j 
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Semiclassical analysis 


Number density & Lee-Yang zeros 
Lee-Yang zeros of Za: 


e o% 
. è 
r . Re 
. e 


Semiclassical study gives the correct transition! 


" 2 (1<p/U « 3/2), 
Nel = 39, Za + 1 (0 « p/U <1), 
Pap 0 (-1/2<p/U <0). 


(YT, Hidaka, Hayata, arXiv:1509.07146[hep-th]) 
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Semiclassical analysis 


Important interference among multiple thimbles 


Let us consider a “phase-quenched” multi-thimble approximation: 


Za) =D le | = e 9 45(0, e 0727). 


o Lee-Yang zeros cannot appear at y/U = 0, 1. 
e One-thimble, or “phase-quenched”, result: n ~ u/U + 1/2. 


Consequence 


In order to describe the step functions, we need interference of 
complex phases among different Lefschetz thimbles. 


(cf. Particle Productions: Dumulu, Dunne, PRL 104 250402) 
(cf. Hidden Topological Angles: Behtash, Sulejmanpasic, Schafer, Unsal, PRL 115 041601) 
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Semiclassical analysis 


Numerical results 


Results for GU = 30: (1, 3, 5-thimble approx.: Jo, Jo U J: 


p1, and Jo U J41 U Jae ) 
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(YT, Hidaka, Hayata, arXiv:1509.07146[hep-th]) 
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Semiclassical analysis 


Numerical results 
Results for GU = 100: 


n 
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Necessary numbers of Lefschetz thimbles ox 5U/27. 
(YT, Hidaka, Hayata, arXiv:1509.07146[hep-th]) 
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Summary & Conclusion 


Summary and conclusion 
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Summary & Conclusion 


Summary and Conclusion 


e Picard—Lefschetz theory gives a suitable framework for 
saddle-point analysis even if S(@) takes complex values. 


e One-site Hubbard model is a nice toy model to play with the 
sign problem. 


o Destructive interference among complex phases of Lefschetz 
thimbles play a pivotal role for the baryon Silver Blaze. 
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